Observation results by the TAMA300 detector 
on gravitational wave bursts from stellar-core collapses 
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We present data-analysis schemes and results of observations with the TAMA300 gravitational- 
wave detector, targeting burst signals from stellar-core collapse events. In analyses for burst grav- 
itational waves, the detection and fake-reduction schemes are different from well-investigated ones 
for a chirp-wave analysis, because precise waveform templates are not available. We used an excess- 
power filter for the extraction of gravitational-wave candidates, and developed two methods for the 
reduction of fake events caused by non-stationary noises of the detector. These analysis schemes 
were applied to real data from the TAMA300 interferometric gravitational wave detector. As a 
result, fake events were reduced by a factor of about 1000 in the best cases. The resultant event 
candidates were interpreted from an astronomical viewpoint. We set an upper limit of 2.2 x 10^ 
events/sec on the burst gravitational- wave event rate in our Galaxy with a confidence level of 90%. 
This work sets a milestone and prospects on the search for burst gravitational waves, by establishing 
an analysis scheme for the observation data from an interferometric gravitational wave detector. 

PACS numbers: 04.80.Nn, 07.05.Kf, 95.85.Sz, 95.55.Ym 



I. INTRODUCTION 

Direct observations of gravitational waves (GWs) will 
reveal new aspects of the universe Since GWs are 
emitted by the bulk motion of matter, and are hardly ab- 
sorbed or scattered, they have a potential to carry astro- 
physical and cosmological information different from that 
by electromagnetic waves. In order to create a new field 
of GW astronomy, several groups around the world are 
developing and operating GW detectors. Among them, 
much effort is being made recently on interferometric de- 
tectors: LIGO 2] in U.S. A VIRGO [3 and GEO 4] 
in Europe, and TAMA Q in Japan. These detec- 
tors have wide-frequency-band sensitivity between about 
10 Hz and a few kHz range, and have an ability to observe 
the waveform of a GW, which would contain astrophysi- 
cal information. In these detectors, both high sensitivity 
and high stability are required because GW signals are 
expected to be extremely faint and rare. 

There are several kinds of target GW sources in 
these interferometric detectors 0, , and data-analysis 
schemes are being developed and applied to the observa- 
tion data. Since GW signals are considered to be faint, 
an efficient data-analysis scheme is required to extract 
the remains of GWs from noisy detector outputs. The 
target GWs are classified by the signal waveforms: chirp 
waves, continuous waves, burst waves, and so on. A chirp 
wave is a sinusoidal waveform with increasing frequency 
and amplitude in time, which is radiated from an inspi- 
raling compact binary just before its merger. Since this 
waveform is well-predicted, an effective and sophisticated 
method of matched filtering is used in chirp- wave search; 
correlations between the data and a template (the pre- 
dicted waveform) are calculated to extract a GW signal 
from noisy data ^ 0| . A continuous wave has a sinu- 
soidal waveform with a stable frequency for over many 
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years, which is radiated from a quasi-stationary com- 
pact binary or a rotating neutron star. A matched fil- 
tering method can also be used in the search for contin- 
uous waves from known sources, because we can predict 
their waveforms precisely by observations with electro- 
magnetic waves prlll^. 

A burst wave, which is the target of this article, is also 
a promising gravitational-waveform class. This wave has 
a short spike-like waveform with a duration time of less 
than 100 msec, which is emitted from stellar-core collapse 
in a supernova explosion or a merger of a binary system. 
Unlike the chirp and continuous-wave cases, a matched 
filtering method cannot be used in a burst-wave analy- 
sis. This is because a set of precise waveform templates 
that cover the source parameters is not available, while 
typical waveforms are obtained by numerical simulations 
El El 113 Thus, several signal-extraction methods, 
called 'burst filters', have been proposed for the detec- 
tion of these burst gravitational waves: an excess power 
filter j3| , a cluster filter in the time-frequency plane ^3 , 
a slope (or a finear-fit) filter [13, and a pulse correlation 
filter ha. Since we have only a little knowledge on the 
waveforms, these filters look for unusual events in the 
Gaussian-noise background. 

For the detection of GWs, evaluation and reduction 
of fake-event backgrounds are critical problems. In each 
analysis scheme described above, we define an evalua- 
tion filter (such as a correlation with the template in a 
matched filtering method, and certain statistics describ- 
ing any unusual behavior of the data in burst analyses) , 
and record the filter output as a GW event trigger if it 
is above a given threshold. The event triggers usually 
contain fake events, which are caused by statistical and 
externally induced excesses in the detector noise level. 
Although an ideal interferometric detector would have a 
stationary Gaussian- noise behavior, the detector output 
is far from stationary in practice, affected by external dis- 
turbances, such as seismic motion, acoustic disturbance, 
changes in the temperature and pressure and so on. As a 
result, it is likely that real signals are buried in these fake 
events, or are dismissed with a larger detection threshold 
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set to reduce fakes. Thus, rejection of these fake events, 
or a veto in other words, is indispensable for the detec- 
tion of GWs. In a chirp-wave analysis case, the output 
of the matched filter is less affected by non-stationary 
noises because it is only sensitive to a waveform similar 
to GWs. In addition, since we know a precise waveform 
of the target GW signal, we can reject the fake events 
by evaluating how well the candidate waveform fits to 
the template 0, 0, . On the other hand, the affects 
of fake events are more serious and critical in the burst 
analysis case. Since burst filters are designed to extract 
any unusual behavior of the detector output, they are 
also sensitive to non-stationary noises by their nature. 
Moreover, it is not straightforward to distinguish these 
fakes from a real signal, and to reject them, because we 
do not know the precise GW waveforms. 

There are several schemes to reject these fake events: 
coincidences by multiple detectors, veto analyses with de- 
tector monitor signals, rejection by waveform behaviors, 
and so on. Among them, the most powerful and reliable 
way will be a coincidence analysis with multiple inde- 
pendent detectors 21, 22, 23J. If we detect gravitational- 
wave candidates with multiple detectors simultaneously 
(or within an acceptable time difference), we can declare 
the detection of a real signal with high confidence. In 
a rough estimation, the fake rate is reduced by a power 
of the number of independent detectors. On the other 
hand, fake reduction with a single detector is also im- 
portant, even in a coincidence analysis as the rejection of 
fakes with a single detector would reduce accidental coin- 
cidences. In observation runs, many auxiliary signals are 
recorded together with the GW signal-channel in order 
to monitor the detector status. Since some of them are 
sensitive to detector instabilities, it is possible to reject 
non-stationary noises with them US- In addition, 
even without precise GW waveforms, fake events are re- 
jected by investigations of the signal behavior with our 
knowledge or assumptions on the waveforms j26| . 

In this article, we present a data-analysis scheme for 
burst GWs, and results obtained by applying them to 
real observation data. The data used in this work were 
over 2700 hours of data obtained during the sixth, eighth, 
and ninth data-taking runs (DT6, DT8, and DT9, respec- 
tively) of the TAMA300 detector Ulllll. We adopted 
an excess power filter as a burst filter, which is robust 
for uncertainties of the GW waveforms Tg', '2^ . In ad- 
dition, we used two fake-reduction methods. One was a 
veto with detector monitor signals. Another was our new 
method of rejection based on the waveform behavior of 
the time scale. Although there havc^been several previ- 
ous works on similar veto methods |2^|^|2^, they were 
applied to a limited subset of observation data. We im- 
plemented these veto methods in our burst analysis code, 
analyzed real observation data, and evaluated their effec- 
tiveness. Such a full-scale analysis is important because 
the effectiveness of the vetoes strongly depends on the 
quality of the real data. 

The obtained event triggers were interpreted from an 
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FIG. 1: Waveform examples from the DFM gravitational 
waveform catalog. The amplitudes are for an event at lOkpc 
distance from the detector. 



astronomical point of view; the results were used to set 
upper limits on Galactic events. We carried out Monte- 
Carlo simulations of Galactic events with waveforms ob- 
tained by numerical simulations of stellar-core collapses. 
In previous works, upper limits by real observations have 
been set with artificial waveform models of short spikes, 
Gaussian waves, or sine- Gaussian waves |2l|,|2^l23- On 
the other hand, realistic waveforms by numerical simula- 
tions have been used to evaluate the efficiencies of burst 
filters with simulated Gaussian noises 0, 0| . We in- 
tended to set upper limits in a realistic way: using a 
realistic distribution of the Galactic events, targeting at 
realistic waveforms obtained by numerical simulations, 
and analyzing long observation data from the detector. 
As a result, we expect to obtain prospects for both cur- 
rent and future research. 

This article is organized as follows. In Section II, we 
overview our burst analysis: target waves, analyzed data, 
and our burst filter. In Section III, we describe veto meth- 
ods with a monitor signal and the signal behavior. After 
that, we present analysis results and an interpretation of 
the results, setting an upper limit on the Galactic events, 
in Section IV. At last, we present discussions and a con- 
clusion of our research in Section V and VI. 



II. GENERATION OF EVENT TRIGGERS 

In this Section, we overview the target GW signal char- 
acteristics, used data, and a burst filter used in this work. 



A. Target gravitational waves 

The target of the analysis in this work is a burst GW 
from stellar-core collapse (core-collapse supernova explo- 
sion). It is difficult to predict its waveform analytically, 
because of the complex time evolution of the mass den- 
sities in the explosion process. Thus, the explosion pro- 
cess and radiated GWs have been investigated by numer- 
ical simulations 0, Q, . Although these simulations 
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Frequenqr [Hz] 

FIG. 2: RSS amplitude and center frequency calculated from 
waveforms in the DFM catalog. The amplitudes are for events 
at the Galactic center (closed circles), and for events at 100 pc 
distance from the detector (open circles) . Each error bar indi- 
cates the frequency range within which the spectrum density 
value is above half of its peak value. The noise level of TAMA 
at DT9 and the design sensitivity of LCGT 36] are shown to- 
gether. 



were performed with differently simplified models, simi- 
lar waveforms were obtained in these simulations. 

Among these simulations, Dimmelmeier et al. have 
presented rather systematic surveys on GWs from stellar- 
core collapses 13]. They have obtained 26 waveforms 
with relativistic simulations of rotational supernova core 
collapses, with axisymmetric models with different ini- 
tial conditions in a differential rotation, an initial rota- 
tion rate, and an adiabatic index at subnuclear densities. 
Although the simulation did not cover all of the initial 
parameters, we used them as reference waveforms in our 
analysis, assuming that typical characteristics and behav- 
ior of the GWs from stellar-core collapses are included in 
this waveform catalog. 

We processed the original waveforms of the catalog 
(we call it DFM catalog in this article) with a 30 Hz 
second-order high-pass digital filter, and resampled them 
to 20 kHz in order to be compatible with the data from 
the detector (described in the next part). Figurenshows 
examples from the waveform catalog. While these wave- 
forms have different behaviors, they have common char- 
acteristics; about a 1 mscc-short spike, and a total dura- 
tion of less than 100 msec. According to the DFM cata- 
log, the averaged amplitude of GWs radiated by super- 
novae at the Galactic center (8.5 kpc distance from the 
detector) is (/ipcak) = 1-5 x 10~^° in a peak strain ampli- 
tude, or (h^-Bs) = 4 X 10^^^ [Hz^^/^] in root-sum-square 
(RSS) amplitude. Here, a RSS amplitude is defined by 



-.1/2 



(1) 



1.2 kHz (Fig. 121), which is around the observation band 
of interferometric detectors. Also, it is estimated from 
the DFM catalog that a total energy radiated as GWs in 



one event is (i^tot) = 8 x 10 
Mq is the mass of the Sun. 



[M0C^], in average. Here, 



B. Data from a gravitational wave detector 
TAMA300 



We applied our analysis method to observation data 
obtained by TAMA300 0, Q; TAMA300 is a Japanese 
laser-interferometric gravitational wave detector, lo- 
cated at the Mitaka campus of the National As- 
tronomical Observatory of Japan (NAOJ) in Tokyo 
(35°40'N, 139°32'E). TAMA300 has an optical config- 
uration of a Michelson interferometer with 300 m- length 
Fabry-Perot arm cavities and with power recycling to en- 
hance the laser power in the interferometer. During the 
operation, the mirrors of the detector are shaken by a 
625 Hz sinusoidal signal, which enables us to calibrate 
the detector sensitivity continuously with a relative er- 
ror of less than 1% 31] . The main output signal of the 
detector, which would contain GW signals, is recorded 
with a 20 kHz, 16 bit data-acquisition system 32]. Be- 
sides the main output signal, over 150 monitor signals are 
also recorded during the observation: signals for the laser 
power in and from the interferometer, detector control- 
loop signals, seismic and acoustic monitor signals, sig- 
nals for temperature and pressure monitor, and so on 
[27I ]. These monitor signals are used for diagnosing the 
detector condition, and for veto analyses (Section 
The recorded data are stored in DLT tapes on site, and 
are sent to data-analysis computers at the collaborating 
institutes by Giga-bit optical network connections. 

In TAMA, nine observation runs have been carried out 
so far since the first observation run in 1999, and over 
2700 hours of data have been collected. In this work, we 
used the data in the sixth, eighth, and ninth data-taking 
runs (DT6, DT8, and DT9, respectively, TablePl. We ob- 
tained over 1000 hours of data in each of DT6 and DT8, 
operating the detector stably and with a good duty cycle. 
The duty cycle in DT9 was not very high, because most of 
the day time was spent for the adjustment and measure- 
ment of the detector during the first half term. On the 
other hand, we obtained data with uniform quality with 



TABLE I: Summary of long data-taking runs by TAMA300. 
The noise level and total observation data amount are de- 
scribed. The last column (D.C.) represents the duty cycle 
throughout the data-taking run. 

Term Noise level Total data D.C. 

[Hz-i/^] [hours] 



where h{t) is the strain amplitude of the GW [23 • The 
central frequencies of the waves range from 90 Hz to 



DT6 Aug. - Sept, 2001 5 x IQ-^ 
DT8 Feb. - April, 2003 3 x 10" 
DT9 Nov., 2003 - Jan., 2004 2 x 10" 



1038 
1157 
558 



87% 
81% 
54% 
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Frequency [Hz] 



FIG. 3: Typical noise level during data-taking runs with the 
TAMA300 detector. The noise level has been improved run 
by run. The spectrum contains several line peaks: harmonics 
of 50 Hz AC power line, violin mode peaks of the suspension 
wire of the mirror (described as 'V'), and a calibration peak 
at 625 Hz (described as 'C'). 
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a high duty cycle in the second half, which included that 
during the quiet time of new-year hoHdays. The duty 
cycle was 96% (207 hours' observation in 216 hours) in 
this quiet term of DT9. We used only this stable term as 
the DT9 data in the following analyses. 

Typical noise spectra in these observation runs are 
shown in Fig.|21 The noise level has gradually been im- 
proved by detector investigations between these runs. 
The detector has a floor-level sensitivity of around from 
300 Hz to 2 kHz frequency range. The floor level is 
2 X 10-21 [Hz-i/2] in DT9 at around IkHz. The spec- 
trum contains several line peaks: harmonics of a 50 Hz 
AC line, violin mode peaks (around 520 Hz and integer 
multiples) of the suspension wire of the mirror, and a cal- 
ibration peak. Since these lines could affect the analysis 
results, they were removed in the data analyses. 



C. Extraction of signals by an excess-power filter 

We developed and implemented a burst-wave analy- 
sis code based on an excess-power burst filter. Am ong 
several filters proposed so far, an excess-power filter [13 
and a TF-cluster filter look for an increase of power 
in the data of a detector, while a slope filter and a 
pulse correlation filter 'l9| monitor correlations between 
the data and assumed waveforms. Roughly speaking, a 
higher detection efficiency is attained with assumptions 
on the waveform. On the contrary, the efficiency is dras- 
tically degraded if there are any errors in the assumed 
waveforms. An excess-power filter is robust because it 
uses only a little information on the target waveforms: 
the signal duration time and the frequency band. The 
evaluation parameter is the total noise power in a given 
time-frequency region. In spite of its robustness, it is 
nearly as efficient as matched filtering for signals with 
short duration and a limited frequency band jlq . 

Our filter generates event triggers in the following steps 



FIG. 4: Data-processing chart of our excess-power burst filter. 
We first calculate a spectrogram from the detector output 
data. Next, we obtain the time-series SNR by averaging the 
frequency components. Then, we extract event triggers by a 
given threshold. 



(Fig.^ details are described in Appendix IXjl : (i) A spec- 
trogram (time-domain change in noise power spectrum) 
is calculated from the output data of the detector; the 
power spectrum is calculated with a At — 12.8 [msec] 
data chunk using a fast Fourier transform (FFT), which is 
repeated with 1.6 msec time-delays. We used a Hanning 
window in each FFT process to obtain smooth spectra 
and time-series results. Here, each spectrum has a fre- 
quency resolution of 1/ At = 78 [Hz]. Since the original 
data contains many line peaks (AC line peaks in every 
50 Hz, etc.), this low-resolution spectrum is contaminated 
by these peaks. Thus, these peaks are removed from the 
original time-series data before calculating the spectro- 
gram, (ii) In each spectrum, power in pre-selected fre- 
quency bands are averaged so as to obtain a time-series 
of averaged power, P„. Each spectrum is normalized 
(whitened) by the typical noise spectrum within 30min 
before a calculation of the average in the frequency com- 
ponents. As a result, P„ represents the signal-to- noise 
ratio (SNR): the ratio of the averaged signal power to the 
typical noise power in pre-selected time- frequency region. 
In this work, we selected a fixed band of A/ = 2270 [Hz] 
from 230 Hz to 2.5 kHz. (iii) Event triggers are extracted 
if the averaged power is larger than a given threshold, 
Pn > ^th- If unusual signals in the detector output are 
sufficiently large, they will be observed in the filter out- 
put. Continuous excesses above the threshold are clus- 
tered to be a single event. Each event trigger is recorded 
with its parameters: the peak signal power Pcv, the time 
of the event tcv, the duration time above the threshold, 
and so on. 
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Injected Amplitude ( 



FIG. 5; Relation of injected signal amplitude and SNR ob- 
tained by the injection test for the DT9 data. The asterisk 
points and the curve represent the signal-injection results and 
the fitting result, respectively. The points at the lower right 
side of the plot are caused by waveforms in which the signal 
power are concentrated at a low-frequency band; the sensitiv- 
ity of TAMA is worse in these frequency bands (see Fig.|5J . 



data-taking run were fitted by 

1 + {Cut X ' ^rss )^ (2) 

where Cdtk represent the averaged efficiency coefficients 
{x: the data-taking number). From the injection results, 
we obtained the coefficient values: Cdt6 = 2.2 x 10^^, 
Cdt8 = 3.3 X IQi^ and Cdt9 = 8.7 x lO^^. 

An inverse of this coefficient corresponds to the GW 
amplitude with which the signal power is the same as 
the noise power by our filter. Thus, it is interpreted 
as the noise-equivalent amplitude of the GW signal. 
The noise-equivalent GW amplitude was ftrss, noise = 
1.1 x 10^^° [Hz^^/'^] for DT9 with our excess-power fil- 
ter. From the estimation that the averaged GW ampli- 
tude was (/irss) = 4 X 10~^^ [Hz~^/^] for a 8.5 kpc event, 
TAMA has the ability to detect events within around 
300 pc away from Earth with this noise-equivalent am- 
plitude. 



The parameters of the filter, length of the time chunk 
(At) for each FFT, and analysis frequency band (A/) 
were selected to be effective for the reference burst GW 
signals. According to the DFM catalog, the signals have 
short spike-like waveforms, i.e. a short duration and a 
wide frequency band. Although the selected parameters 
(At = 12.8 [msec]. A/ = 2270 [Hz]) were not fully op- 
timized for the waveforms, the analysis results were not 
changed very much with a different parameter set. More- 
over, we could keep the robustness of the excess-power 
filter with this rough tuning of the time-frequency bands. 



D. Signal-injection simulations 

The output of the filter (Pcv) is a dimensionless signal- 
to-noise ratio (SNR). SNR is cahbrated to a physical 
value, such as the GW amplitude, by results of signal- 
injection simulations (called 'software injection tests') 
with real data from the TAMA300 detector. In the sim- 
ulation, target reference waveforms are superimposed to 
the detector data with proper calibration (estimated from 
the transfer functions of the detector, the whitening fil- 
ter, the anti-aliasing filter before the data- acquisition sys- 
tem). The signals are injected and analyzed one by one 
with equal time separations in the data so as to evaluate 
the data uniformly. The amplitudes and waveforms were 
selected randomly from 10^^^ < h^ss < 10"^^ [Hz^^/^j 
and from 26 waveforms in the DFM catalog, respectively. 
This data were analyzed by the same code as that for the 
raw data analysis. 

The results of the signal-injection test are shown in 
Fig-El the recorded power SNRs of the events (Pev) are 
plotted as a function of the root-sum-square amplitudes 
of the injected signal (/irss)- The injection results of each 



III. REDUCTION OF FAKE EVENTS 

In this Section, we describe veto methods to reject fake 
events caused by detector instabilities. We have used two 
veto methods: a veto method using auxiliary signals for 
the detector monitor, and a veto method by the waveform 
behavior: the time-scale of the signal. 

A. Veto with auxiliary signals for the detector 
monitor 

Here, we describe a veto method using auxiliary signals 
recorded together with the main output of the detector: 
a correlation between the monitor signal and the main 
output, the confirmation not to reject real GW signals, 
and a false dismissal rate estimation. 



1. Veto with the intensity monitor signal 

We investigated some of the monitor signals, and found 
strong correlations between the short spikes in the main 
output and the monitor signal for the laser intensity 
in a power-recycling cavity of the interferometer. Fig- 
urelHl shows the correlation between the main output and 
the intensity monitor signal. The intensity-monitor data 
were processed by the same excess-power filter to detect 
short-spike instabilities. The filter parameters were the 
same as that for the main signal analysis, except for the 
frequency range. In the analysis of the intensity-monitor 
data, the frequency range was A/int = 1170 Hz from the 
DC frequency, which was determined from the spectrum 
shape of the burst spikes in the intensity data. The closed 
circles in Fig. represent the power (SNR) of the events 
in the GW data-channel and in the intensity monitor. 
In this figure, the event triggers were extracted with a 
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FIG. 6: Correlation between the powers in the main output 
channel and the intensity monitor channel. The closed cir- 
cles are event candidates selected by the threshold Pth = 4; 
the gray asterisks are powers at the time 100 msec-shifted 
from the event candidates, and the open circles are results 
of hardware-injection tests. 



FIG. 7: Estimation of an accidental coincidence. The 
intensity- veto threshold was decided so that the accidental co- 
incidence probability (or the false-dismissal rate in this case) 
would be 1%. As a result, the thresholds are Pth, int ~ 3.9, 
2.2, and 3.0 for DT6, DT8, and DT9, respectively. 



threshold of Pth ~ 4 for the DT8 data. For a compar- 
ison, the powers at the time 100 msec-shifted from the 
triggers are also plotted in this figure (gray asterisks) 
pnl . This figure shows that there are strong correlations 
between these two signal powers for some of the event 
triggers, and only weak correlations outside of the event 
triggers. 

Thus, vetoes of fake events with the intensity monitor 
signal are expected to be effective; when the outputs of 
two excess power filters (one for the GW signal-channel, 
and another for the intensity monitor) exceed the respec- 
tive thresholds simultaneously, the triggers are labeled as 
fakes, and are removed from the event candidate list. 



2. Estimation of a false-dismissal rate 

To use the veto with the intensity monitor signal, wc 
should confirm that the intensity instabilities were not 
caused by huge GW signals; otherwise, we may reject 
real GW signals by this veto. During DT8, we investi- 
gated the response of the detector by injecting simulated 
waveforms to it. In this test, which is called a 'hardware 
injection' test, we shook the interferometer mirrors with 
a short-spike waveform and a typical burst-waveform ob- 
tained by numerical simulations |l4j with various am- 
plitudes. The results of the hardware-injection test are 
plotted as open circles in Fig.El There were 147 in- 
jections, and 117 events were above the event-selection 
threshold of Pth = 4. We observed no excess power in 
the intensity monitor with an intensity veto threshold of 
-Pth, int — 2.2 (described below), while the injected sig- 
nals appeared with sufficiently large powers in the filter 
output for the GW signal-channel. The number of acci- 
dental excesses with this threshold is expected to be 1.2 
(1% of injected events). Thus, the result that no excess 
was found above the intensity threshold is well-consistent 



with the expected accidental background. From these re- 
sults, we ensured that the intensity instabilities were not 
caused by huge GWs, and that it is safe to use the inten- 
sity monitor signal for a veto analysis. 

The threshold for the intensity excess power is selected 
so as to reduce fake events effectively with small proba- 
bility to reject real GW signals, i.e. with a small false- 
dismissal rate. The false-dismissal rate is equal to the 
accidental coincidence rate between the intensity excess 
and the excess in GW signal-channel, which was esti- 
mated from the distribution of the power in the intensity 
monitor signal. We selected a threshold so that the false- 
dismissal rate would be 1%, which resulted in thresholds 
of Pth,int = 3.9, 2.2, and 3.0 for DT6, DT8, and DT9, 
respectively (Fig.O. The difference in the thresholds for 
these data-taking runs were caused by the difference in 
the detector stability in each run and the improvement of 
the intensity monitor instrument between DT8 and DT9. 

B. Veto by a signal behavior test 

Here, we describe the second veto method, a veto 
method by the time-scale of the signal. Statistics for the 
signal evaluation, and an estimation of the false-dismissal 
rate are described. 



1. Signal behavior test 

As described above, fake events are rejected by care- 
ful selection and an investigation of the monitor signal- 
channels. However, it is hard to see any clear correla- 
tions for all of the fake events in practice, because there 
are various origins of the fakes, which are difficult to be 
identified. Thus, a test of the signal behavior at the 
main output of the detector will be helpful to reduce 
fake events. 



Q 




FIG. 8: Evaluation of a data point on a C2-C1 correlation plot. 
The curve shows the expected locus for a reference waveform, 
which was drawn by sweeping the signal amplitude. The dis- 
tance between a data point and a reference point on the refer- 
ence curve represents the similarity of the signal amplitudes 
and their time scales. 



The effectiveness of the veto with a signal behavior test 
depends on how well we know, or how many assumptions 
we set, on the signal behavior. In the burst-wave analy- 
sis case, the waveforms by numerical simulations suggest 
that GWs from stellar-core collapse have a short dura- 
tion, typically less than 100 msec. We know that some of 
the detector instabilities last longer than a few seconds 
from experience. Thus, some of the fakes caused by these 
slow instabilities are rejected by evaluating the time scale 
of the event triggers |2a |. 

2. Evaluation statistics for the time-scale veto 

In this work, we selected to evaluate the time scale 
of event triggers with statistics around the event (Ap- 
pendixlBl for details). The excess-power filter outputs 
P„ , a time series of the power in a pre-selected frequency 
band. We calculate two statistics from the ±AT/2 time- 
series data around the event-candidate time tav- 

ci = Qi and C2 = ^ - 1^ , (3) 

where Qk = ^Pn/NJk: integer) is the fc-th-order 
moment of the power [s^, [s^. Here, N is the num- 
ber of power-data points in the evaluation time (within 
icv ± Ar/2). The statistics ci, which is related to an av- 
eraged power, has information about the stability of the 
noise level on a given time scale. On the other hand, C2 is 
related to the second-order moment of the noise power. 
Since it is normalized by the averaged power (Qi), the C2 
value becomes constant if the signal power is much larger 
than the background noise level. The constant number 
is determined by the time scale of the event: large in a 
short-burst case, and small in a case of a slow change in 
the noise power. 




C2 value 

FIG. 9: Correlation plot of ci and C2 used for the rejection of 
fake events, calculated with a time window of AT = 0.82 [sec]. 
The solid curves show the expected loci for the reference wave- 
forms from the DFM catalog. Among them, the curve for 
the waveform with the longest time-scale is shown as a bold 
curve. The dotted curves represent the contours for the dis- 
tance Dmin to the longest time-scale reference curve. When 
event triggers are in the gray area, which is the -Dmin > 5 
area, they are rejected as fakes. The closed circles represent 
the results of the signal-injection test using the data of DT9. 
With a fake-event-selection threshold of -Dth = 5, the false 
dismissal rate was estimated to be 0.3%. 



We use an evaluation parameter (i^min)7 which rep- 
resents the similarity to the GW signal, estimated from 
the ci and C2 statistics. When we plot the data point in 
the C2-C1 plane (Fig.|Sl) , they will be in the left region for 
long-duration event cases and in the right region for short 
burst-like event cases, and will be in the upper region for 
large power events. Thus, two independent information 
on each event, the power and the time scale, are shown 
by the position of each data point in this plot. Here, 
the distance {D) between the data point and a reference 
point (expected position of the real GW signal, estimated 
from the signal waveform), represents the similarity of 
their signal behaviors. The evaluation parameter (I?min) 
is the smallest distance for all of the reference waveform 
and amplitude combinations. 

In a practical application of this method to real data, 
we set a loose selection criteria for fake events (Fig.EJ. 
We have 26 reference GW waveforms from the DFM cat- 
alog, which have different time scales. Instead of com- 
paring the time scale of an event trigger with that of 
each reference waveform, we compare it only with the 
longest time scale of the reference GW waveforms. In 
other words, the evaluation parameter of an event (I?rnin) 
is set to be zero if its time scale is shorter than that of 
the longest reference waveform. 

3. Selection of parameters for the time-scale veto 

Fake events with different behaviors from that of the 
real ones are rejected when the minimum distance (Z^min) 
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FIG. 10: Event-trigger rate results with TAMA data. The 
horizontal and vertical axes represent the threshold (Pth) and 
the trigger rate in a unit of events/sec, respectively. The 
analysis result with simulated Gaussian noise is also plotted 
together with the DT6, DT8, and DT9 results. 
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FIG. 11: GW amplitude and corresponding trigger rate; the 
event rate (vertical axis) with larger amplitude than a given 
hrss (horizontal axis) is plotted. 



is larger than a given threshold (-Dth)- Here, two param- 
eters should be set in this veto analysis: an evaluation 
time- window (AT) and a fake-selection threshold (Dth). 
The time window is selected to be AT = 0.82 sec, so 
that the veto would be effective. Since this time win- 
dow length is between the typical time scales of the fakes 
(larger than a few seconds) and the real GW signals (less 
than 100 msec), we can expect a clear distinction between 
them. On the other hand, the threshold for the rejection 
(-Dth) is selected so that the false-dismissal rate of real 
GW signals would be acceptable. 

The false-dismissal rate was directly evaluated from 
the results of a signal-injection test with the real data 
from the TAMA300 detector (described in Sec. InD)! . 
This simulation is important because the real data from 
a detector do not have an ideal Gaussian noise distri- 
bution. The closed circles in Fig.El show the results of 
the signal-injection test to the DT9 data. The injection 
results are distributed well-around the theoretical predic- 
tions shown as solid curves. With a fake-event-selection 
threshold of Dth = 5, the false-dismissal rate was esti- 
mated to be 0.3%: 6 injection results were rejected out 
of 2006 injections MT|. Although this value is larger than 
that estimated by a statistical analysis, it is sufficiently 
small for a veto analysis. The false dismissal rates were 
also investigated for the DT6 and DT8 data with a simi- 
lar signal-injection test, and found to be 0.6% and 2.9%, 
respectively. The differences come from the original be- 
havior of the data in these observation terms. 



IV. DATA-PROCESSING RESULTS WITH THE 
TAMA300 DATA 

The analysis method described above was applied to 
real data from TAMA300. In this section, we consider 
the results of the TAMA data analysis with the vetoes, 
and the interpretation of the results from an astronomical 



point of view. The data were mainly processed by a 
PC-cluster computer placed at the University of Tokyo. 
This machine is comprised of 10 nodes, and has 20 CPUs 
(Athlon MP 2000-t- by AMD Inc.). The analysis time for 
the excess power filter was about 30-times faster than the 
real time; it took about 1 /30 sec to process 1-sec data. 

In the data processing, the first 9-min and the last 1- 
min data of the each continuous observation span were 
not used because they sometimes contained loud noises 
caused by detector instabilities, or excited violin-mode 
fiuctuations. In addition, the duration time of rejected 
fake events is considered as a dead time of the detector, 
and subtracted from the total observation times. The 
dead times by the fake rejections were 1.3%, 1.7%, and 
0.4% of the observation time for DT6, DT8, and DT9, 
respectively. The effective observation times (Tobs) are 
shown in TablellTl 



A. Event-trigger rates 

FiguretTUI shows the event-trigger rates obtained by 
the TAMA data analyses; the trigger rate (in a unit 
of events/sec) is plotted as a function of the event- 
extraction threshold (Pth)- The analysis result with sim- 
ulated Gaussian noise is also plotted in Fig. EH together 
with the DT6, DT8, and DT9 results. Assuming that 
the real GW signals are rare and faint, we can regard 
most of the triggers as being fakes. From this figure, 
one can see that the trigger rates were reduced in these 
data-taking runs with vetoes. For a given GW power 
threshold (-Pth), the trigger rates were reduced by 1/10- 
1/1000. The power threshold could be reduced (for better 
GW-detection efficiency) by a factor of 10-100 for a given 
trigger rate. FiguretTTl shows the event-trigger rates plot- 
ted as a function of ft-rss amplitude, which was obtained 
from Eq. ^ . The detector was gradually improved dur- 
ing the intervals of these data-taking runs. The event- 
trigger rates were reduced from DT6 to DT9 by about a 
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FIG. 12: Reduction factor by two veto methods for DT6, 
DTS, and DT9. The bold curves are the reduction ratio with 
two vetoes. The dotted and dashed curves are only with an 
intensity veto and a time-scale veto, respectively. 



effective than the intensity veto. This is because DT6 
data contained noisy data originating in an instability 
of the laser source and seismic disturbances during the 
daytime. 

The rates are still much larger than that with simu- 
lated Gaussian noises, even with the vetoes and the im- 
provements in the detector. In addition, the trigger rate 
is still much higher than the expected rate of supernova 
explosions. The expected GW event rate is one event 
in a few tens of years, i.e. about 10~^ events/sec in our 
Galaxy. (Here, note that TAMA has an ability to de- 
tect only events within 300 pc away from Earth at best.) 
These results suggest that most of the observed trigger- 
events were fake events caused by an instability of the 
detector, even with vetoes. 



few orders for a given GW amplitude, and by about an 
order for given trigger rates fTablejll)) . 

FiguretT^ shows the reduction factor of event triggers 
with two veto methods; the ratio of the number of event 
triggers after and before the vetoes are plotted as a func- 
tion of the SNR threshold. In these runs, both of the two 
vetoes contributed to the reduction of the rates. They 
worked in complementary ways. The intensity veto is ef- 
fective to short-duration fakes and large SNR fakes. On 
the other hand, the time-scale veto is effective for the 
long-time instability of the detector output and small 
SNR fakes. 

In DT9, many event triggers were rejected as fakes by 
the intensity veto. This is because the pre-amplifier and 
whitening filter for the data acquisition of the intensity 
signal were improved in this run. Since the reduction fac- 
tor is better for large SNR fakes, it is expected that the 
reduction ratio will be further improved with a higher 
detection efficiency of the intensity instabilities. On the 
other hand, only a small fraction of fakes were rejected 
by the time-scale selection in DT9 because the detector 
operation was sufficiently stable. The detector was oper- 
ated very stably thanks to a quiet seismic environment 
during the holiday weeks in the second half of DT9. In 
addition, the drift of the typical noise level was small at 
that time. In DT6, the time-scale veto was much more 



TABLE II: Summary of data analysis results. The noise- 
equivalent GW RSS-amplitudes (^rss, noise), the dead times by 
the vetoes (Trcj), the total effective observation times (Tobs), 
the trigger rates for h^ss > 1 x 10"" [Hz"^/^], and the GW 
RSS-amplitudes above which the trigger rates are one event 
per hour are described. 



/irss, noise Trej Tobs Rate 1-hour amp. 

[Hz-^/^] [hours] [hours] [sec~^] [Hz"^/^] 

DT6 4.5 X 10"^" TO 937.8 2.1 x 10"^ 2.1 x 10"'** 

DTS 3.0 X 10"^° 18.0 1064.2 7.0 x 10"* 1.4 x 10"" 

DT9 1.1 X 10"^° 0.8 194.6 2.5 x 10"^ 2.5 x 10"^'-* 



B. Simulations of Galactic events 

Since the event-trigger rate is still much larger than 
that expected from ideal Gaussian noise or observed su- 
pernova rates, we cannot claim the detection of GW sig- 
nals from the data-analysis results. Thus, we set up- 
per limits for stellar-core collapse events in our Galaxy. 
We carried out Monte-Carlo simulations with a source- 
distribution model of our Galaxy, and with waveforms 
from the DFM catalog. The simulated data were ana- 
lyzed in the same way as the detector data, and compared 
with the observation results. 

In the simulation, we adopted a source-distribution 
model based on the observed luminous star distribution 
in our Galaxy, assuming that the event distribution of 
the stellar-core collapses was identical to it. There have 
been studies on the star-distribution model based on sky- 
survey observations [33 . l34j . In our simulation, we used a 
simple axisymmetric distribution model (an exponential 
disk model) described in a cylindrical coordinates, 

p(i?,0,^)cxexp(^-A_M^ , (4) 

where p, Rq — 3.5 kpc, and ho — 320 pc are the density 
of the events, and the characteristic radius and height of 
the density of the Galactic disk, respectively. As well as 
the non-axisymmetric components, such as spiral arms, 
the thick disk and halo structures were neglected in our 
simulations because their number of stars was only about 
3% of that of the disk component |3^ . We adopted Rq — 
8.5 kpc and — 20 pc for the position of the Sun in our 
simulation. 

We used 200 hours stable data in the second half of 
DT9 for the Galactic signal-injection test. This test was 
performed according to the following steps: (i) Set the 
GPS times at which simulated events are injected; these 
times are uniformly separated between the start and end 
times of the observation run. Decide the position of each 
event randomly according to the Galactic-event distri- 
bution described by Eq. I^J . Select the waveform of each 
event randomly from the DFM catalog, (ii) Calculate the 



11 



>. 
o 
c 
.9! 
o 

ie 
a 
c 
o 



10 



10"- 



£ 10"' 



10"' 

















- 


ie.l 

B 

3 1 
3 




'•••...^ Event 
**'•'.,. (a 
''*'•• 


ate in DT9 


O, 

S' 


Effici 


/\ 

ency for 


**» 


' Threshold 





10" 



10"^ g 
m 

0) 

Hn-2 C 
10 0) 
> 

10"^ aj 

(0 



10"' 



10 10 10 

Event Power Threshold {P,h) 



10'' 



FIG. 13: Results of a Galactic-event simulation. The solid 
curve shows the detection efficiency for the Galactic events 
as a function of the threshold. The trigger rate obtained by 
the DT9 observation is also plotted as a dotted curve for a 
comparison (right vertical axis). 



distance and sky direction seen from the detector for each 
event, from the position of the event in the Galaxy and 
the injection time information, (iii) The expected GW 
amplitude is calculated from the distance to the event 
source and the detector antenna pattern for the sky po- 
sition of the event We assumed non-polarized GWs; 
the GW power is equally distributed to the two polariza- 
tions, (iv) Inject each event waveform to the TAMA300 
data with estimated amplitude, and analyze the data 
with a similar code as that for the raw-data analysis, 
(v) Extract the events at the injected time. 



C. Results of Galactic-event simulations 

FiguretTSl shows the results of a Galactic-event simu- 
lation. The fraction of the detectable Galactic events 
(the detection efficiency, egai, left axis) is plotted as a 
function of the SNR threshold (Pth). The event-trigger 
rate in DT9 is also plotted for a comparison (right 
axis). With an event-selection threshold of Pth = 3.0 
(which corresponds to averaged amplitude of /irss.th = 
1.6 X 10-20 [Hz-i/2] for DT9), the detection efficiency 
was estimated to be Cgai = 3.1 x 10^^ for the Galactic 
events. The threshold was selected so that the expected 
contribution of the Gaussian noise would be sufficiently 
small (less than 1% of the triggers above the threshold). 
The upper limit for the event rate determined from the 
TAMA raw-data was -Rdt9,ul = 6.8 x 10"^ events/sec 
with a confidence level of 90%. From these results, 
we obtained the upper limit for the Galactic event rate 
to be i?gai,uL = i?DT9,UL/egai = 2.2 X 10^ [events/sec] 
with a 90% confidence level This value is consid- 

erably larger than the theoretical expectation of about 
10"^ events/sec. 

Besides the upper limit for the rate of a stellar-core 
collapse in our Galaxy, an upper limit was set for the GW 



energy rate. The total energy radiated as GW, i^tot, was 
estimated for each event from its waveform. The upper 
limit for the energy rate radiated as GW was estimated 
by the product of the event-rate upper limit, i?gai,uL, 
and the averaged GW energy of the events, {Etot)- As 
a result, we obtained ii^cw.UL — 3.8 x 10^* [Mqc^ /sec]. 
Again, this value is considerably large; the rate of the 
total energy radiated as GWs would be about Mq^x/ (2 x 
lO'^) [Mqc^ /years], where Mcai is the total mass of our 
Galaxy, which we assume to be 2 x 10^^ M©. 

There are uncertainties in setting the upper limits by 
several origins. Here, we consider the effect of the de- 
tector calibration error, statistical error in the Monte- 
Carlo simulation for the Galactic events, and the error 
in the Galactic model. The calibration error in the con- 
version of the detector output to the GW strain ampli- 
tude was estimated to be less than 1%. This calibration 
error causes an amplitude error in the signal-injection 
test. From Eq. (|2Jl and the results of the Galactic signal- 
injection test, the uncertainty in the upper limit is es- 
timated to be 2.9% with a detector calibration error of 
1%. On the other hand, the statistical error in the Monte- 
Carlo simulation is determined by the number of the sim- 
ulated events above the threshold. We generated 3.8 x 10® 
Galactic events, and detected 1.2 x 10"* events above the 
threshold of Pth — 3.0. Assuming that the number of 
detected events follows a Poisson distribution, the event- 
rate uncertainty is 0.9%. At last, the error in the Galac- 
tic model would affect the results. The parameter i?o in 
the Galactic model has an error of 9.4% |3^]. This er- 
ror corresponds to a 5.8% uncertainty in the detection 
efficiency for the Galactic events (cgai) with a threshold 
of Pth = 3.0, which was estimated by additional simula- 
tions. In total, the uncertainty in our upper-limit results 
is 6.5% at most. The detection efficiency will be reduced 
(the upper limit will be increased) by 1.2% by including 
a thick disk and halo components. 



V. DISCUSSIONS 
A. Comparison with previous studies 

We interpreted the observation results from the view- 
point of the Galactic event rate in the previous section. 
In this part, we interpret the results in an similar way as 
in the previous studies for a comparison; we set an upper 
limit on the rate of GW events incident on the detector 
as a function of the GW amplitude '2^, "SO^ . 

We used 200 hours of data in the DT9 stable term, and 
set an upper limit on the rate of the events received by the 
detector, following the procedure to set the upper limits 
by the LIGO group 29J. At first, the detection threshold 
was fixed, and the upper limit was set at the number of 
the events above the threshold with a given confidence 
level. Then, the detection efficiency for a given GW am- 
plitude was estimated by a software signal-injection test. 
Here, the events were distributed randomly on the celes- 
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FIG. 14: Detection efficiency and an upper limit (90% con- 
fidence level) for the waveforms of the DFM catalog. The 
detection efficiency with a given threshold of Pth = 10* is 
shown in the upper plot, and the upper limit on the event 
rate is shown in the lower plot (solid curve). The gray dotted 
curve is the envelope of the upper limits with various detec- 
tion thresholds. 



FIG. 15: Detection efficiencies and upper limits (90% con- 
fidence level) for the Gaussian (time scale of 0.5 and 
1 msec) and sine-Gaussian (central frequency of 554, 850, and 
1304 Hz) waveforms. The gray dotted curve is the envelope 
of the upper limits for the 1304 Hz sine-Gaussian waveform 
with various detection thresholds. 



B. Outlook for the detection of burst GWs 



tial sphere in order to include the directivity of the de- 
tector. We investigated the Gaussian (with a time scale 
of T — 0.5 and 1 [msec]) and sine-Gaussian waveforms 
(with Q-value of 9 and central frequencies of 554, 850, 
and 1304 Hz) for a comparison with the previous studies, 
as well as the waveforms from the DFM catalog. Finally, 
we estimated the upper limit on the event rate from the 
upper limit on the number of events, A'^ul, the detec- 
tion efficiencies, e, and the observation time, Tobs, by 
i?UL = A^UL/(e-Tobs). 

We set the threshold to be Pth = 10^, which resulted 
in one trigger above the threshold. Assuming Poisson 
statistics, we obtained the corresponding upper limit of 
3.89 events with a confidence level of 90%. The detec- 
tion efficiencies and the upper limit results are shown in 
Fig-El (waveforms from the DFM catalog) and Fig.ll5l 
(Gaussian and sine-Gaussian waveforms). 

The upper limit for sufficiently large events (ex. hrss > 
1 X 10~^^ [Hz^^/^]) was 0.49 events/day with a confi- 
dence level of 90%. This upper limit is comparable to the 
LIGO-Sl result of 1.6 event/day [2^ and the Glasgow- 
MPQ coincidence result of 0.89 events/day [23|. On the 
other hand, the resonant detector network has set an 
upper limit of ~ 4 x 10"'^ events/day j30|. These differ- 
ences in the upper limits come mostly from the differ- 
ent observation times. As for the sensitivity for smaller 
amplitude signals, the GW amplitude for 50% detection 
efficiencies (averaged over source directions) was around 
1 X 10~^^ Hz~^/^ for short and high-frequency waveforms 
in our case. The upper limit curve is almost comparable 
with the LIGO-Sl results for high-frequency signals, and 
larger for lower frequency (< 800Hz) ones p9l |. 



The large event rate and upper limit results show that 
the detector output is still dominated by fake events, even 
after these vetoes. Thus, further research efforts are nec- 
essary to detect burst gravitational waves. In this part, 
we discuss the outlook to better vetoes, coincidence anal- 
yses with other observatories, and better performance of 
the detector and data-processing scheme. 

In Section inra we presented a veto analysis method 
with only one monitor signal, an intensity monitor. Sim- 
ilar methods can be used with the other monitor signals 
along with careful investigations of their correlations with 
the main output of the detector. However, we have found 
no other monitor signal with a clear correlation so far: 
laser power at the signal-port (the dark port), monitors 
for the seismic fluctuations, an acoustic monitor signal. 
Thus, it is necessary to investigate more deeply the mon- 
itor signals, and to introduce better monitor signals that 
are sensitive to the detector instabilities. 

There are other event-selection criteria than the time- 
scale selection method presented in Section llll Bl For ex- 
ample, the time scale of an event can be simply evaluated 
by the duration time above the event-selection threshold. 
In this case, we should consider that the veto results will 
be strongly dependent on the event amplitude ^4^. We 
will be able to reduce fake events even further by knowing 
the common characteristics of the target events, and set- 
ting them as event-selection criteria. For this, more sys- 
tematic and precise simulations of stellar-core collapses 
and investigation on the waveform will be helpful. 

Coincidence analyses with the other detectors for 
GWs, electromagnetic waves, and neutrinos will improve 
the result significantly, though we have focused on the re- 
duction of fakes with a single detector in this article. The 
observation runs by TAMA300 (DT8 and DT9) were car- 
ried out at the same term as the LIGO second and third 
scientific (observation) runs (called S2 and S3), and co- 



incidence analyses are underway |35|- We note that our 
work discussed in this article is also a part of the LIGO- 
TAMA coincidence analysis; the list of the event triggers 
obtained in our work will be used in the coincidence anal- 
ysis. 

In addition to a reduction of fakes, the improvements 
of the detector both in the floor noise level and in the 
reduction of non-stationary noises are also important. 
The performance of the TAMA300 detector has grad- 
ually been improved from DT6 to DT9 concerning both 
the noise level and the stability, and the detector still 
has room for improvement. In addition, burst filters with 
higher efficiencies are under development in the TAMA 
group and other groups. Since we can only observe events 
within about the 300 pc range with the current sensitiv- 
ity of TAMA, the detection efficiency for the Galactic 
events is very small (cgai = 8.9 x 10~^ with a threshold 
for a noise-equivalent GW amplitude). The sensitivity 
should be improved by about two orders so as to cover 
our Galaxy, and to realize a sufficiently large detection 
efficiency. This sensitivity will be realized by the next- 
generation detectors, such as LCGT (Fig.|3) ,3^ and ad- 
vanced LIGO 113. 



VI. CONCLUSION 

We presented data-analysis schemes and results of ob- 
servation data by TAMA300, targeting at burst signals 
from stellar-core collapses. Since precise waveforms are 
not available for burst gravitational waves, the detec- 
tion schemes (the construction of a detection filter and 
the rejection of fake events) are different from those for 
chirp wave analyses. We investigated two methods for 
the reduction of non-stationary noises, and applied them 
to real data from the TAMA300 interferometric gravi- 
tational wave detector. As a result, these veto meth- 
ods, a veto with a detector monitor signal and a veto by 
time-scale selection, worked efficiently in a complimen- 
tary way. The former and the latter were effective for 
short-spike noises and for slow instabilities of the detec- 
tor, respectively. The fake-event rate was reduced by a 
factor of about 1000 in the best case. 

The obtained event-trigger rate was interpreted from 
the viewpoint of the burst gravitational-wave events in 
our Galaxy. From the observation and analysis results, 
we set an upper limit for the Galactic event rate to be 
2.2 X lO^^ events/sec (confidence level 90%), based on 
a Galactic disk model '34l and waveforms obtained by 
numerical simulations of stellar-core collapses [T5 |. In 
addition, we determined the upper limit for the rate of 
the energy radiated as gravitational-wave bursts to be 
3.8 X lO-^Mgc^/sec (confidence level 90%). These large 
upper limits show that the detector output was still dom- 
inated by fake events, even after the selection of events, 
and gives us prospects on both current and future re- 
search: the necessity for further improvement of the anal- 
ysis schemes, coincidence analyses with multiple detec- 
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FIG. 16: Example of the line removal results for the DT9 data. 
Time-series data before (plotted in gray) and after (plotted in 
black) the line removal are shown. The lower plot is a zoom 
up of the spike in the upper plot. 



tors, better predictions on the waveforms, and future de- 
tectors, such as LCGT and advanced LIGO, to cover the 
whole of our Galaxy. This work has set, we believe, a 
milestone for these research activities. 
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APPENDIX A: DETAILED DESCRIPTION ON 
THE EXCESS-POWER FILTER 

1. An excess-power filter 

In this part, we detail the excess-power filter. We as- 
sume that the output of a detector is comprised of an 
ideal stationary Gaussian noise, nit), and a signal, s{t) 
(non-Gaussian component caused by gravitational waves 
or instability of the detector): v{t) = n{t) + s{t). The 
power spectrum is calculated for every At data chunk 
with given time-delays (St), using a fast Fourier trans- 
form (FFT). As a result, we obtain a spectrogram (a 
time-frequency plane) of the noise (or signal) power with 
dt time separation and 1/ At frequency resolution. Here, 
the Fourier component is also described by the sum of the 
noise and the signal: Vmn = fimn + Smn, where m and n 
represent the indices for the frequency and time, respec- 
tively. Then, the power in each time-frequency compo- 
nent is described by \vmn\^ ■ 

In order to make our filter equally effective for all of the 
analysis frequencies, we normalize the power spectrum by 
an averaged noise level, 7V„ = Yl°=no''"~^ |nmnP/A^av, 
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where A'av is the number of time components used for 
the average. Then, the normahzed power is written as 

where Nmn = |"-mnP/_/V_m, Smn = l^mnP/ Njn, and 

Cmn = ^{smn ' "-mn}/ ^m, meaning a normahzed noise 
power, a normahzed signal power, and a normahzed cor- 
relation between the signal and noise components, re- 
spectively. Since fimn has a Gaussian noise distribution, 
Njnn has a distribution of two degrees of freedom (an 
exponential distribution): P{N„in) = eKp(—Nmn)- 

Then, the output of the excess-power filter, the aver- 
aged power for a given time component, is written as 




Frequency [Hz] 



FIG. 17: Results of line removal for the DT9 data. The noise 
spectrum before (plotted in gray) and after (plotted in black) 
the line removal are shown. 



m 

where M is the number of frequency components used in 
the average; only the values of the power in pre-selected 
frequency components are used to calculate the averaged 
power, Pn- From Eq. IjAip . P„ is written as 

Pn = Nn + 2Cn + Sn, (A3) 

where = N^n/M, C„ ee C„„/Af, and Sn = 
J2m Smn/M. This represents the time evolution of the 
power in detector output in given frequency bands. 

2. Data conditioning 

The data from the detector is not an ideal Gaussian 
noise, in practice; the noise spectrum is not white, the 
noise level changes in time, and many line peaks are in- 
cluded in the output. Thus, data conditioning before 
processing the excess power filter is indispensable. 

The spectrum contains several line peaks: harmonics 
of 50 Hz AC line, violin mode peaks (around 520 Hz and 
integer multiples) of the suspension wire of the mirror, 
and a calibration peak. These line peaks are removed 
in the following processes. At first, we obtain a Fourier 
spectrum from 72 sec of data by FFT. We then set the 
line-frequency components to be zero. In addition, the 
lower frequency components below 160 Hz are also re- 
jected. At last, we obtain a time-series data by calculat- 
ing the inverse FFT of the spectrum. With this process, 
the line peaks are clearly removed from the spectrum 
(Fig.^1 I17|l . Moreover, since only a small number of 
frequency components are rejected, the burst waveforms 
are not vary much affected by the line removal process. 

The frequency and time dependences of the noise spec- 
trum are compensated by the normalization with the av- 
eraged noise power spectrum, Nm- We calculated by 
averaging the power spectra for 30 min before the data 
analyzed by the excess-power filter. In order to avoid 
the large spikes from disturbing the averaged spectrum, 
we rejected noisy 0.7% spectra from those used in each 



average. We found that we could obtain stable aver- 
aged spectra, and that each spectrum was normalized 
well with this method. 

APPENDIX B: TIME-SCALE EVALUATION AND 
VETO 

1. Evaluation in AT time chunk 

In this part, we describe the details of the veto method 
with time-scale evaluation of the event triggers. In our 
veto methods, each event is evaluated by the statistics 
in a tev ± AT/2 data chunk {tcv- the time of the event). 
Here, data points of the excess-powcr-filter output are 
contained in the time window AT, i.e. AT = NSt. From 
the output of the excess-power filter, P„, we define the 
evaluation parameters ci and C2 as 

Ci =Ql, C2 = y^- 1, (Bl) 

where Qi, and Q2 are the first- and second-order mo- 
ments of Pn for N data points, respectively, written as: 

^ no + N-l no+W-l 

Ql^- Y: Q^^^ E (Pr.r. (B2) 

n— no n—UQ 

Here, note that Qi is an averaged power for M x time- 
frequency components. On the other hand, C2 is defined 
by the second-order moment normalized by the averaged 
power. This value is analogous to the kurtosis (defined by 
the fourth-order moment of data), which describes any 
non-Gaussianity of the data [s^ |33 . 

2. Statistics of Qi and Q2 

We calculate the statistics of parameters Qi and Q2 as 
a preparation for calculating the statistics of the param- 
eters ci and C2 , defined in Eq. (|B1|I . Although the M xN 
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components are not independent in practice, because of 
overlapping in time and a window function, we describe 
the following calculations while assuming that they are 
independent, for simplicity. In practical use, the statis- 
tics (the averages, the variances, and the covariance) are 
estimated by replacing the time and frequency window 
size, N and M, by effective time and frequency range, 
A'eff and Moff. The effective window sizes, Ncs and Mcs, 
are estimated by simulations with Gaussian noises. 

From Eq. HA3|I . Qi is also written by the sum of the 
noise, signal, and their correlation terms. Thus, the ex- 
pected value is 



EiQi) 



E{Nn) + 2E{Cn)} = a + l, (B3) 



where we define the average of the signal component 
power by 



a ■ 



(B4) 



and we use relations E{Nn) = 1 and E{Cn) = 0. The 
expected value of the square of Qi is written as 

ml) = eI^Y^Y. (^J- + ^0 + 2Q) {Si + Ni+ 2Ci, 



2a + I 
MN 



+ (a + l)^ 



where we use relations E{NkNi) = E{Nk)E{Ni) [k ^ I), 
E{Nl) = (M + 1)/Af , and so on. Thus, the variance of 
Qi is written as 



On the other hand, the expected value of Q2 is 
E{Q2) = -B |^l^(5„ + iV„ + 2C„)'^ 



(B5) 



= /32Q; +2a + l 



2a -f 1 



(B6) 



where [32 is a constant value related to the second-order 
moment of the signal. 



(320" ^Y^n/N. 



(B7) 



Similarly, a constant value, (3^, is written as S'^/N = 
/Saa^. The constant numbers a, P2 and are deter- 
mined only by the waveform and the amplitude of the 
signal. The value a represents the normalized signal 
power. The value P2 depends on the time scale of the 
signal; P2 becomes large for a short signal. 

With more complicated, but similar, calculations, the 
variance of Q2 and the covariance between Qi and Q2 



are obtained to be 



Mii(Qi:Q2) 



8 ^ o 20M + 32^ 2 

16(M2 + M + 3) 
^ M37V " 
2(2M2 -t- 5M -f 3) 

^ APN ' 

2/32a^ 



MN 



3(Af + 1) M 
^ -a • 



M M 



Statistics for Ci and C2 



(B8) 



From the results described above, we obtain the ex- 
pected value and the variance of ci as 



E{ci) = a + l, 
2a + 1 



M2(ci) 



MN 



(B9) 



On the other hand, the expected value and the variance 
, of C2 are obtained as 

E{C2) = H, + o(J^ 

fi2ic2) = ^i2{Ql)Hf + 2/in(Qi, Q2)HiH2 + fi2{Q2)Hl 
1 



+0 



(BIO) 



where 



Ha = C2{EiQi),EiQ2)), 
H2 = 



dQi 
dc2 



Qi=E(Qt),Q2 = E(Q2) 



dQ2 



Qi=E{Qi),Q2 = E{Q2). 



Thus, we obtain the mean and variance of C2, 
{^2 - l)a^ + ^ 



E{C2) 
M2(C2) 



(a + 1)2 
4a2 



(Bll) 



MN{a + 1)6 
{2(/33 - I3l)a^ + (4/33 - ^Pi - ^2)^^ 
+2(/33-/32)a+(/32-l)}, 

and the covariance of c\ and C2, 

2(/32 - \)o? 



(B12) 



Mii(ci, C2) 



MN{a 1)3 



(B13) 



Here, we have neglected the higher terms, such as 
0(1/M27V), 0[\IM^N). 
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FIG. 18: Theoretical predictions of the data point on a C2- 
ci plot for given waveform parameters {a and /32)- The 
loci corresponds to the /32 parameter of /32 ~ 122.21 x 
(0.2,0.4,0.6,0.8,1). The locus for a waveform with a large 
/32 value (with a short time scale) appears on the right side 
of the plot. 



4. Veto method 

From the calculations described above, we can esti- 
mate the expected values as E{ci) and E{c2) when the 
waveform and amplitude are given. Figure lTHl shows the 
expected points for given waveforms in the C2-C1 plane; 
each curve is plotted by sweeping the power a. When a 
is small, the (c2, ci) point is around (0, 1) independently 
of the waveform {/32)- (Here, we assume M ^ 1.) On 
the other hand, the position has a strong dependence on 
the /32 parameter for large a: C2 — > /32 — 1, for a ^ 1. 
Since signals with different time scales have different (32 
values, they appear along different loci. 

Since the average time, AT, is finite, and since the 
data contains the Gaussian noise components, the (c2, 
ci) data points have a distribution around the predicted 
curve shown in Fig. 1181 With an approximation that AT 
(or data point number N) is sufficiently large, the (c2, 
ci) points for given a and P2 have a two-dimensional 
Gaussian distribution, determined by the expected val- 



ues {E{ci) and E{c2)), the variances (/X2(ci) and /i2(c2)), 
and the covariance ^ii(ci,C2). We define the distance 
between a reference point, which is calculated by a grav- 
itational waveform (a reference waveform), and a data 
point by 

= ^{fl2{c2)Acl 

-2/iii(ci, C2)AciAc2 -I- Ai2(ci)Ac^} , (B14) 
where Aci, Ac2, and V are defined by Aci = ci — i?(ci), 

AC2 = C2 - £'(02), and V = ^i2{Cl)^J,2{C2) - Aill(ci,C2), 

respectively (Fig.lSJ. This distance, which is normalized 
by the variances and covariance, represents how similar 
these points are; if D is small, the data point has a sim- 
ilar amplitude and a time scale as that of the reference 
point. Here, the minimum distance (-Dmin) for various 
amplitudes (a) is the similarity of the time scales of the 
data point and the reference waveform. Thus, we use 
-Dmin to distinguish fakes from the true GW signals. If 
the data points have a two-dimensional Gaussian distri- 
bution around their mean point, and if the variances and 
covariance are sufficiently smaller than the curvature of 
the locus, the minimal distance (Z^min) approximately 
obeys an exponential distribution, 

PiD^in) « e-^»."/^ (B15) 

For a practical implementation of the veto scheme, we 
should consider the acceptability of multiple reference 
waveforms and a reduction of the computational load. 
Thus, we adopted a conservative way for the veto analy- 
sis; we reject only the events with longer time scales than 
the longest one in the reference waveforms (Fig.O. We 
use one waveforms with the smallest P2 value, i.e. with 
the longest time scale, as the reference waveform, and set 
-Dmin = if the data point is on the right side of the ref- 
erence locus in the C1-C2 plot. In addition, we prepare a 
map of Dmin in the C1-C2 plane before the data analysis, 
so as to reduce the computational load during the data 
analysis. 
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